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ABSTRACT 



^ [ We present the discovery of a radio galaxy at a likely redshift of z = 4.424 in one 

of the flanking fields of the Hubble Deep Field. Radio observations with the VLA and 
MERLIN centered on the HDF yielded a complete sample of microjansky radio sources, 
of which about 20% have no optical counterpart to I < 25 mag. In this Letter, we 
address the possible nature of one of these sources, through deep Hubble Space Telescope 
NICMOS images in the FllOW (Jno) and F160W (Hieo) filters. VLA J123642+621331 
has a single emission line at 6595 A, which we identify with Lya at z = 4.424. We argue 
that this faint (Higo = 23.9 mag), compact {r^ ^ 0'.'2), red (Igu — K = 2.0) object is 
most likely a dusty, star- forming galaxy with an embedded active nucleus. 

Subject headings: galaxies: active — galaxies: starburst — galaxies: evolution — galax- 
ies: individual (VLA J123642+621331) 



^Based on observations with the NASA/ESA Hubble Space Telescope obtained at the Space Telescope Science 
Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract 
NAS5-26555; and with the W. M. Keck Observatory, which is operated as a scientific partnership among the University 
of California, the California Institute of Technology, and NASA, made possible by the generous financial support of 
the W. M. Keck Foundation. 



1. Introduction 

One of the contemporary topics of interest in extragalactic astronomy is to measure the star 
formation history of the universe (Madau et al. 1996). However, a major uncertainty in calculating 
the star formation rate (SFR) is the effect of dust obscuration, particularly for measurements in the 
ultraviolet and at high redshifts. For example, Calzetti & Heckman (1999) found that at z ~ 3 even 
a modest amount of dust can reduce the 1500 A flux of a galaxy by a factor of 5. This is of particular 
significance for the high-z evolution of the SFR: does the SFR turn over at z ~ 1.5 as Madau et al. 
(1996) suggested, or is it constant (or even increasing) for z > 1.5 (e.g. Pascarelle, Lanzetta, &: 
Fernandez-Soto 1998; Steidel et al. 1999)? If a high-redshift starburst is significantly reddened by 
dust, then its ultraviolet flux will be so obscured that it will have a negligible Lyman- limit break. 
Such sources would not have been detected by Steidel et al. (1999), and their measurements of the 
SFR could still be under-estimating the actual amount of star formation at high-z. 

A powerful technique to avoid the problem of dust obscuration is to use radio-selected galaxies 
as probes of the star formation history (Cram 1998). The most important advantage of this 
technique is that radio emission is not attenuated by dust and thus a radio-selected sample is 
unbiased with respect to dust, unlike the optical/infrared selection techniques that are widely used. 
Haarsma et al. (1999) conclude from their study that the radio-selected SFR is somewhat higher 
than the Madau et al. (1996) results for z < 1, even when the latter have been corrected for dust. 
At higher redshifts, the "radio Madau diagram" is poorly constrained but is consistent with a 
relatively constant or increasing star formation rate. This suggests that a significant fraction of the 
star formation in the universe may be obscured by dust. 

Radio-selected samples thus provide an essential tool for understanding the cosmic history of 
dust, as well as of the SFR. In particular, radio sources that are heavily reddened in the optical 
(restframe ultraviolet at redshifts of interest) can be used to constrain the dust content of high-z 
galaxies. We present in this Letter deep HST/NICMOS images of the microjansky radio source 
VLA J123642-I-621331. At a probable redshift of 4.424, this is the second highest redshift radio- 
selected galaxy currently known. It was barely detected in the reddest optical wavebands and is 
an excellent candidate for a dusty high-z galaxy. 

We use a cosmology with Hq = 65 km s^-*^ Mpc^-^, Qm = 0.2 and J7a = throughout. All 
magnitudes are in the AB system, unless noted otherwise. For comparison, Vega-based magnitudes 
are approximately given by: I8i4 — 0.4, Jno — 0.7, Higg — 1-3, K — 1.9, and (I814 — K) -|- 1.5. 



2. VLA J123642+621331 

Two of the deepest high-frequency radio surveys available are those centered on the Hubble 
Deep Field (Richards et al. 1998; Richards 2000) and the Small Selected Area 13 (Windhorst et al. 
1995). These surveys reach 5a detection hmits of 8 //Jy at 8.5 GHz (HDF and SSA13) and 40 /nJy 



at 1.4 GHz (HDF). Optical identifications of the sources were made on deep Hubble Space Telescope 
and ground-based images, primarily in the I-band. Approximately 20% of the 1.4 GHz sources were 
not identified, although a few of them show faint optical emission below the formal completeness 
limit of the I-band images. The full sample is discussed by Richards et al. (1999); in this Letter we 
investigate in detail one of these "unidentified" sources. 

VLA J123642-I-621331 is a steep-spectrum radio source (a = 0.94 it 0.06) with flux densities of 
70 /iJy at 8.5 GHz and 470 //Jy at 1.4 GHz (Richards et al. 1998; Richards 2000). Approximately 
10% of its radio flux density is extended in an eastward jet in OdS MERLIN observations, with the 
remainder in a resolved compact core (Muxlow et al. 1999). The source was identified with a faint 
object below the formal completeness limit of the HDF Flanking Field I-band image (IW3), with 
a magnitude of I8i4 = 25.3 =b 0.2 (Figure la). The object was also identified in the K-band image 
of Dickinson (1998) from the KPNO 4-meter telescope, and has a magnitude of K = 23.25 it 0.05 
(Figure Id). The KPNO images do not detect the object in either the J or H bands. Aussel 
et al. (1999) identify VLA J123642-I-621331 in their supplementary list of sources detected by 
ISOCAM on the Infrared Space Observatory. It has a flux of 23J]]^2 ^Jy ^^ the LW3 filter at 15 ^m, 
corresponding to an AB magnitude of 20.5 ± 0.5. The source lies at the edge of the SCUBA sub- 
millimeter map of Hughes et al. (1998) and was not detected, giving a conservative upper limit of 
5 mJy to its 850 ;um flux. 



3. NICMOS Observations & Processing 

In December 1997, we observed VLA J123642-I-621331 with the Hubble Space Telescope NIC- 
MOS camera 2 for three orbits in FllOW (close to the J-band) and six in F160W (essentially the 
H-band), in the Continuous Viewing Zone. In each orbit, five 1024-second exposures were taken 
using a spiral dither pattern with l'.'3 offsets. We reduced the images using a two stage process: 
first we ran them through the standard pipeline CALNICA, using the most recent calibration files, 
and then we applied a further flat-field correction to the data, adapted from standard ground-based 
infrared imaging methods. 

Following CALNICA processing, additional bad pixels were flagged after a visual inspection 
of the data — these consisted of a dead column, the coronograph hole and other insensitive pixels 
(also known as the "grot"). Since the NICMOS camera actually consists of four physically separate 
sub-arrays, we obtained better results by dividing each partially-reduced image into four separate 
quadrants. For each filter, all the exposures of each quadrant were stacked and a median image 
calculated. Given that the exposures were dithered by 1'.'3, this produced a map of the residual 
instrumental features, devoid of any astronomical objects — a "super-sky". The super-sky images 
for each quadrant were normalized to the mean of all four quadrants in order to preserve the 
quadrant-to-quadrant photometric accuracy. Each exposure was then divided by this normalized 
super-sky. We also tested the results of subtracting a scaled copy of the super-sky from each 
exposure and found that the results were not significantly different (the standard deviation of 



an empty sky region varied randomly by less than 3% between the two methods). The success 
of this super-sky division suggests that the instrumental features left after CALNICA processing 
were most likely due to differences between the sky and the calibration flat-fields, rather than poor 
subtraction of the bias or dark current. Finally, all the quadrants were combined to produce FllOW 
and F160W mosaics. By ensuring that the mean of our super-sky was unity, the overall sensitivity 
of the images was unchanged by the additional processing and they were calibrated using the most 
recent photometric calibration parameters from STScI. 

Figure Ic shows the 6-orbit NICMOS image in F160W. The counterpart to the radio source is 
clearly detected at 1.6 /um, with an AB magnitude of Higo = 23.87 it 0.04 in a 1.5 arcsec diameter 
circular aperture. The 3-orbit NICMOS image in FllOW provided a marginal detection of Jno = 
25.2±0.4 (Figure lb). Thus the object is very red with (Ism-K) = 2.0±0.2, (Jno-Hieo) = 1.3±0.4 
and (Hieo-K) = 0.6 ±0.1. 



4. Keck Observations 

We obtained spectra of VLA J123642+621331 through l'.'5 wide, « 30" long slits using the 
LRIS spectrograph (Oke et al. 1995) at the Keck II telescope in slitmask mode. On UT 1998 
February 19, with 0'.'8 seeing and photometric conditions, we observed VLA J123642+621331 for 
1.9 hr (position angle 103°) with the 400 lines mm^^ grating, sampling the wavelength range 5700- 
9400 A, at AAfwhm ~ H A. On UT 1999 May 10, with 0'.'6 seeing and thin cirrus, we observed 
the source for 2 hr (position angle —67°) with the 150 lines mm~-^ grating, with AApwHM ~ 17 A 
resolution over the wavelength range 4000 A to 1 /xm. Between each 1800 s exposure, we performed 
a ~ 4" spatial shift along the slit to facilitate removal of fringing in the reddest regions of the 
spectra. Final wavelength calibration is accurate to better than 1 A. 

There is a strong, single emission line at A ~ 6595 A in both data sets, which we identify with 
Lya at a redshift of 4.424 (Figure 2, Table 1). The line shifts by ~ 7 A between the two observations, 
which corresponds to a velocity difference of 320 km s^^ if it is due to a Doppler shift. In both years' 
data, the line was offset by ~ 1" to the north-west of the (marginal) I-band detection. Emission-line 
regions of high-redshift radio galaxies are known to be kinematically complex (Chambers, Miley, & 
van Breugel 1990; van Ojik et al. 1997), thus slight pointing changes between the two observations 
may have caused the slit to sample different regions of spatially-extended, line-emitting gas. 



5. Discussion 

The only two reasonable identifications for the single emission line are Lya or [O II] 3727. We 
argue that it is unlikely to be [O II] at z = 0.77 for the following reasons. The faint K magnitude 
of 23.3 argues strongly against the source being at low-z when it is compared to the K-z relation 
of radio galaxies (see van Breugel et al. 1999 for a recent version). If VLA J123642+621331 were 



at z = 0.77, it would be 3-4 magnitudes under luminous in K compared with all other known 
radio galaxies at that redshift. The restframe equivalent width if the line were [O II] would be 
WrQjji > 207 A, which is very large for [O II]. The absence of a redshifted [O III] doublet at 
~8860 A further argues against z = 0.77, although this argument is not completely satisfying as 
galaxies show a wide range in [O III]/[0 II] ratios and our flux limits are not particularly strong. 

The alternative is to identify the emission line with Lya at z = 4.424. We compared the 
observed SED from 0.8 //m to 15 ^m with the 1996 revision of the spectral evolution models of 
Bruzual & Chariot (1993). We used a single burst model of solar metalicity and added a foreground 
screen of dust, modeling the effects of dust obscuration with the extinction law of Calzetti, Kinney, 
& Storchi-Bergmann (1994). The best-fitting model we obtain has an age of (1.6 it 1.0) x 10^ years 
with A{V) = 1.6 lb 0.3 mag (Figure 3). No model of any age is able to reproduce the red colors 
of this galaxy with A{V) < 0.5 mag at the 99.99% confidence level. The results are essentially 
independent of the metalicity, from 0.02 of solar to solar. We note that the models predict a 15 ;um 
flux that is too faint, although given the large ISO errors the models were consistent to within 
< 2a. We also attempted to fit the galaxy with a model spectrum at z = 0.77 and were unable to 
obtain any acceptable fit to the optical/near-infrared SED: the minimum x^ was ~ 30 for a 18 Gyr 
model with no dust extinction, compared with a x^ of 2.5 for the best-fitting z = 4.424 model. In 
particular, the low redshift models could not reproduce the red Hieo — K color of 0.6 it 0.1, nor did 
they fit the 15 //m flux. 

The sub-millimeter (850 fim) to radio (1.4 GHz) spectral index has recently been used by 
Carilli & Yun (1999) as a redshift indicator. With a spectral index of af^^ > —0.4 (Si, oc z^^"), 
VLA J123642-I-621331 is consistent with this correlation at either z = 0.77 for a starburst galaxy, 
or at ^ = 4.424 if the 1.4 GHz flux is dominated by an AGN. Thus the current SCUBA sub- 
millimeter detection limit does not favor either a high or low redshift interpretation. The lack of 
a far-infrared/sub-millimeter detection prevents us from drawing conclusions about the emission 
properties of the dust that must be present based on the reddening of the optical/infrared SED. 

We note that it is possible that we detected a serendipitous z = 4.42 background object in 
the Keck spectra, that is not associated with the radio source. We used the number counts and 
redshift distribution of galaxies in the HDF to find the random probability that a z > 4 source with 
a magnitude 1-2 mag below the completeness limit of the I-band image would be found within a 
l" radius of the optical identification. From this calculation we estimate that there is a probability 
of < 10^^ that the spectra are of a background source. 

The surface brightness profile of VLA J123642+621331 was produced by fitting elliptical 
isophotes to the F160W image (Figure 4). We extracted model surface brightness profiles from an 
ensemble of two-dimensional exponential and de Vaucouleurs models convolved with a model PSF 
from "Tiny Tim". The best-fitting model is an exponential with scale-length re = 0'.'l7ib0'.'02 (where 
I oc e^*"' ''") or 1.4ib0.2 kpc in our chosen cosmology. Any unresolved point source (i.e., an AGN) con- 
tributes < 10% of the total flux in the H-band. A de Vaucouleurs r^'^ profile does not fit the data. 



We note briefly the contrast between the dominant exponential profile of VLA J123642+621331 
and the profiles of other very red high-z sources, of which several are now known at z ~ 1.5-2 that 
have regular r^'^ profiles (e.g., Stiavelli et al. 1999, Waddington et al. 1999). 

It can be seen in Figure 1 that there is some two-dimensional structure to the galaxy. It is 
extended in a roughly north-south direction in the F160W image, perpendicular to the radio jet. 
In the F814W image there is a faint hint of extended emission to the west of the radio position, 
in the opposite direction to the radio jet but aligned with the Lya emission. The marginal J-band 
detection may be similarly extended to the west. The apparent structure in the K-band is probably 
not real, but is likely due to the PSF and/or the drizzling procedure, given that other objects in 
the image show a similar structure. It is also possible that there is a contribution to the K-band 
light from redshifted [O II] 3727. Most high-z radio galaxies have distorted or aligned morphologies 
(e.g.. Best, Longair, & Rottgering 1997) and it is an important question whether such a low power 
radio source as VLA J123642+621331 is also aligned with its ultraviolet/optical emission. 

The preceding arguments lead us to the following interpretation. VLA J123642+621331 is 
a disk system at a likely redshift of z = 4.424, containing an embedded, weak AGN. With an 
8.4 GHz luminosity of 2.0 x 10^^ W Hz~^, it is at the faint end of the radio luminosity function 
for AGN, intermediate between an AGN and a starburst. Both the radio jet that is visible in the 
MERLIN image and the sub-millimeter to radio spectral index suggest it is an AGN. The AGN 
must be obscured in the ultraviolet/optical, as it is not seen in either the F160W profile or in the 
colors of the SED. The observed Lya luminosity of 2 x 10*^^ erg s~^ (Table 1) is a factor of 10-100 
fainter than typical high-redshift radio galaxies (Rottgering et al. 1997), and its Lya line- width 
(FWHM ~ 440 km s^^) is 2-3 times narrower than the other known z > A radio galaxies (van 
Breugel et al. 1999). 

The galaxy underwent a burst of star formation approximately 16 million years before we 
observe it, although the presence of dust suggests that this was not the first starburst. It is a dusty 
galaxy with an extinction of 1.6 mag in the visual, and a corresponding ultraviolet extinction of 
A(1216 A) ~ 4 mag. This would be sufficient to render any Lya emission from the disk undetectable 
in our Keck spectrum, if it were of comparable intrinsic luminosity to the (unreddened) line that 
we do observe. The Lya emission line came from an extended region about 1" to the north-west of 
the galaxy, roughly aligned with the radio axis. This Lya emission region could be an infalling gas 
cloud that is either scattering/reradiating Lya from the AGN into our line of sight, or is the site of 
recent (jet-induced) star formation. Such extended Lya clouds have been observed around several 
high-redshift AGN at z ~ 2.4 (Francis, Woodgate Sz Danks 1997; Windhorst, Keel, & Pascarelle 
1998; Keel et al. 1999) and we should perhaps expect such structures to be more abundant at higher 
redshifts, where the gas has had less time to collapse. Deep imaging of the source at ~ 6600 A is 
needed in order to investigate the nature of this cloud. 

We thank Eric Richards, Ken Kellerman, Ed Fomalont and Tom Muxlow for useful discussions 
and for sharing unpublished data; Andy Bunker for contributing to the Keck observations; and 



Mark Dickinson for making his KPNO infrared observations of the HDF pubhcly available. This 
work was supported by NASA grant GO-7452.0*.96A from STScI under NASA contract NAS5- 
26555, and by NSF grant AST9802963. 



REFERENCES 

Aussel, H., Cesarsky, C. J., Elbaz, D., k Starck, J. L. 1999, A&A, 342, 313 

Best, P. N., Longair, M. S., & Rottgering, H. J. 1997, MNRAS, 292, 758 

Bruzual A., G., & Chariot S. 1993, ApJ, 405, 538 

Cazetti, D., & Heckman, T. M. 1999, ApJ, 519, 27 

Calzetti, D., Kinney, A. L., & Storchi-Bergmann, T. 1994, ApJ, 429, 582 

Cariin, C. L., & Yun, M. S. 1999, ApJ, 513, L13 

Chambers, K., Miley, C, k van Breugel, W. 1990, ApJ, 363, 21 

Cram, L. 1998, ApJ, 506, L85 

Dickinson, M. 1998, in The Hubble Deep Field, eds. M. Livio, M. Fall, k P. Madau, Cambridge 
Univ. Press, p. 219 

Francis, P. J., Woodgate, B. E., k Danks, A. C. 1997, ApJ, 482, L25 

Haarsma, D. B., et al. 1999, in preparation 

Hughes, D. H., et al. 1998, Nature, 394, 241 

Keel, W. C, Cohen, S. H., Windhorst, R. A., & Waddington, I. 1999, AJ, in press (December 1999) 

Madau, P., Ferguson, H. C, Dickinson, M. E., Giavalisco, M., Steidel, C. C, k Fruchter, A. 1996, 
MNRAS, 283, 1388 

Muxlow, T., et al. 1999, in preparation 

Oke, J. B., et al. 1995, PASP, 107, 375 

Pascarelle, S. M., Lanzetta, K. M., Fernandez-Soto, A. 1998, ApJ, 508, LI 

Richards, E. A. 2000, ApJ, in press 

Richards, E. A., Fomalont, E. B., Kellerman, K. I., Partridge, R. B., Windhorst, R. A., Cowie, L. 
L., k Barger, A. J. 1999, ApJ, this volume 

Richards, E. A., Kellermann, K., Fomalont, E., Windhorst, R., k Partridge, R. 1998, AJ, 116, 103 



Rottgering, H., van Ojik, R., Chambers, K., van Breugel, W., & de Koff, S. 1997, A&A, 326, 505 

Steidel, C. C, Adelberger, K. L., Giavalisco, M., Dickinson, M., k Pettini, M. 1999, ApJ, 519, 1 

Stiavelli, M., et al. 1999, A&A, 343, L25 

van Breugel, W., de Breuck, C, Stanford, S. A., Stern, D., Rottgering, H., k Miley, G. 1999, ApJ, 
518, L61 

van Ojik, R., Rottgering, H. J. A., Miley, G. K., k Hunstead, R. W. 1997, A&A, 317, 358 

Waddington, I., Windhorst, R. A., Peacock, J. A., Dunlop, J. S., Cohen, S. H., k McClure, R. 
1999, in preparation 

Windhorst, R. A., et al. 1995, Nature, 375, 471 

Windhorst, R. A., Keel, W. C, k Pascarelle, S. M. 1998, ApJ, 494, L27 



This preprint was prepared with the AAS lATpjX macros v5.0. 




Fig. 1. — Multiband images of 

VLA J123642+621331. (a) WFPC2 F814W (I- 
band) grey-scale image, overlaid with the 1.4 GHz 
combined (uniform weighted) MERLIN/VLA radio 
image of Muxlow et al. (1999), at a resolution of 
approximately 0.15 arcsec. Contour levels are at 
flux densities of 8, 16, 32, 128 /iJy (ct = 4 ^Jy). (b) 
The source is barely detected in a 3-orbit NICMOS 
FllOW (J-band) image, (c) The 6-orbit NICMOS 
F160W (H-band) image shows a very red object at 
the radio position, (d) Ground-based K-band image 
with the KPNO 4-meter (I'.'O seeing). Each image 
is 4 arcsec on a side and has been smoothed with a 
2-pixel wide gaussian. 



10 



1150 



Rest Wavelength (A) 
1200 1250 



1300 




6200 6400 6600 8800 
Observed Wavelength (A) 



7000 



Fig. 2.— Keck spectra of VLA J123642+621331. 
The single emission line is detected in two indepen- 
dent observations over a period of more than a year, 
and is identified as Lya at z = 4.424. The 1998 spec- 
trum is smoothed with a 3-pixel boxcar filter. 
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Fig. 3. — Optical/infrared spectral energy distri- 
bution of the radio galaxy. The three spectra are 
the best-fitting model (solid line), the oldest & least 
dusty model (dot-dash line) and youngest & most 
dusty model (dotted line), both at the 3-0" confidence 
limits. 
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Fig. 4.— The NICMOS F160W surface brightness 
profile of VLA J123642+621331. Triangles denote 1- 
a upper limits. The best-fitting model (solid line) is 
an exponential with scale-length 0'.'17 (1.4 kpc), con- 
volved with the appropriate PSF (dotted line) from 
"Tiny Tim". 



Table 1. Spectroscopic Measurements 



Parameter 



UT 1998 Feb 19 UT 1999 May 10 



I(A)''''''''''''''''''''''^^^^^ 

Fhya (10 -'^'' erg cm ^ s ■^) 

^cont (;lo-21 erg cm-2 s"i A^^) 

wi^liA) 

FWHMLya (km s^^) 

LLya (10^2 erg s-1) 



4.427 ±0.001 

6597.4 ± 1.0 

0.74 ±0.14 

-0.2 ±7.3 

>483 
420 ± 75 
2.4 ±0.4 



4.421 ± 0.001 
6590.4 ± 1.1 
0.48 ±0.11 

8.4 ±4.8 
> 249 

464 ± 155 

1.5 ±0.4 



